SUMS OF LAPLACE EIGENVALUES ROTATIONS AND 
TIGHT FRAMES IN HIGHER DIMENSIONS 



R. S. LAUGESEN AND B. A. SIUDEJA 

Abstract. The sum of the first n > 1 eigenvalues of the Laplacian is 
shown to be maximal among simplexes for the regular simplex (the reg- 
ular tetrahedron, in three dimensions), maximal among parallelepipeds 
for the hypercube, and maximal among ellipsoids for the ball, provided 
the volume and moment of inertia of an "inverse" body are suitably 
normalized. This result holds for Dirichlet, Robin and Neumann eigen- 
values. Additionally, the cubical torus is shown to be maximal among 
flat tori. 

The method of proof involves tight frames for euclidean space gener- 
ated by the orbits of the rotation group of the extremal domain. 

The ball is conjectured to maximize sums of Neumann eigenvalues 
among general convex domains, provided the volume and moment of 
inertia of the polar dual of the domain are suitably normalized. 



1. Introduction 

In this paper we show how the Method of Rotations and Tight Frames 
developed in our earlier paper for plane domains [20] can be extended to 
yield higher dimensional bounds on sums of eigenvalues of the Laplacian. 
The transition to higher dimensions requires a broader algebraic context, 
in terms of irreducible symmetry groups. Furthermore, the results must be 
reformulated in a non-obvious way, because the obvious generalization fails 
in every dimension higher than two. 

We will obtain sharp bounds on the sum of the first n > 1 eigenval- 
ues of linear images of highly symmetric domains. This work generalizes 
Polya's results [23, 24] in three ways. His proof could handle only the first 
Dirichlet eigenvalue of planar domains; our method bounds eigenvalue sums 
of arbitrary length, in arbitrary dimensions, and under any major bound- 
ary condition (Dirichlet, Neumann or Robin). For a detailed comparison of 
Polya's approach with our own, see the end of Section 6 in [20]. 

Write Ai, A2, ... for the Dirichlet eigenvalues of the Laplacian on a domain 
in M. d . Let V be the volume and I the second moment of mass around 
the centroid (which in two dimensions simply equals the moment of inertia 
about a perpendicular axis through the centroid). We show that for any 
linear transformation T and any domain D possessing sufficient rotational 
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symmetry, the normalized eigenvalue sum 

yl+2/d 



(Al + • • • + \ n )V 2 l d 



T(D) I 



(1.1) 

T~ 1 (D) 



is maximal when T is the identity matrix. That is, the original domain D 
is the maximizer. 

Both factors in expression (1.1) are scale- and rotation- invariant. Hence 
maximality holds also when T is any multiple of an orthogonal matrix. 

Interestingly, the moment of mass normalization in (1.1) is imposed on 
the "inverse domain" T^ 1 (D), rather than on the domain T(D) where the 
eigenvalues are computed. This feature seems unavoidable in dimensions 
three and higher, as the rectangular box example reveals in Section 3. In 
two dimensions one can normalize the moment of inertia on T(D) [20]. 

We obtain similar results for Neumann and Robin boundary conditions, 
and for flat tori. The Robin case requires a significant extension of the proof 
from the planar situation. 

There are at least two good reasons for studying eigenvalue sums. The 
first is physical: the sum of eigenvalues represents the energy needed to 
fill the lowest n quantum states under the Pauli exclusion principle. The 
second is mathematical: summability methods provide a tools for studying 
high eigenvalues, which are difficult to study individually. An example of the 
gains possible by summation is that while Polya's conjecture (claiming the 
Weyl asymptotic is a lower bound for each Dirichlet eigenvalue) is still open, 
Li and Yau [21] have shown that sums of eigenvalues are indeed bounded 
below by the analogous Weyl asymptotic. 

Notice our work in this paper is geometrically sharp, with an extremal 
domain existing for each fixed n. The Li- Yau inequality is not geometrically 
sharp, but is asymptotically sharp — equality holds for each domain in the 
limit as n — > oo. Asymptotically sharp upper bounds in the Neumann 
case were obtained by Kroger [19], and sums of functions of eigenvalues 
were studied recently by Frank, Geisinger, Harrell, Hermi, Laptev, Loss and 
Weidl ([7, 8, 11, 14] and references therein). 

Eigenvalue estimates for rotationally symmetric domains were previously 
studied by Ashbaugh and Benguria [1], while some higher dimensional upper 
bounds for the first Dirichlet eigenvalue were proved by Freitas and Krejcifi'k 
[9]. Their result applies to general domains, however it is usually weaker 
for the domains we can handle. Nevertheless, it is equivalent for the first 
Dirichlet eigenvalue of ellipsoids. Eigenvalues of certain simplexes have been 
found explicitly [16, 18, 25], although the regular tetrahedron is not one of 
them because the trigonometric method of Lame for the equilateral triangle 
fails to extend to higher dimensions. 

For a detailed literature survey regarding eigenvalue estimates, one could 
begin with our two-dimensional paper [20] and Henrot's book [15]. 
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2. Assumptions and notation 

Eigenvalues. For a bounded domain D in R d , d > 2, we denote the Dirich- 
let eigenvalues of the Laplacian by Xj(D), the Robin eigenvalues by Pj{D), 
and the Neumann eigenvalues by pj(D). In the Robin and Neumann cases, 
we assume the domain has Lipschitz boundary so that the spectrum is dis- 
crete. Denoting the eigenfunctions by uj in each case, we have 

— Auj = XjUj in D j — Auj = pjUj in D j — Aiij = PjUj in D 

uj = on dD \^ + cruj = on dD [ ^=0 on9D 

and 

< Ai < A 2 < A 3 < . . . < pi < p 2 < p 3 < • • • = pi < p 2 < P3 < ■ ■ ■ 

Geometric quantities. Let 

V = volume, 

/ = second moment of mass about the centroid. 
That is, 

1(D) = I \x - x\ 2 dx 
Jd 

where the centroid is x = J D x dx/V . 

Note that / can be interpreted as an average moment of inertia with re- 
spect to a randomly chosen axis (see for example [10]). It can also be viewed 
as the moment of inertia of D in with respect to the axis perpendicular 

to the hyperplane containing D and passing through the centroid. 

Given a matrix M, write its Hilbert-Schmidt norm as 

\m\Hs = (E^) 1/2 = ( trMtM ) 1/2 , 

where denotes the transposed matrix. 

Tight frames. The leading role in our proofs is played by the concept of a 
tight frame: 

A set of unit vectors {yi}f = i in M d forms a unit-norm tight frame if for 
every vector z G IR rf , 

This property resembles the Plancherel identity for an orthonormal basis, 
except that the number of vectors N can exceed the dimension d. 

Tight frames have become an important tool in applied harmonic analysis, 
in recent years. The overdetermined representations they provide are useful 
for noise reduction and robustness to erasures, in signal processing. Note 
that the tight frames in this paper consist of equal- norm vectors; for more on 
that special case, see the work of Benedetto and Fickus [3] and Casazza and 
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Kovacevic [4]. General tight frames do not impose the equal- norm restric- 
tion. All tight frames arise as (rescaled) projections of orthonormal bases 
in higher dimensional spaces [13, Chapter 5]. For this and more informa- 
tion about frame theory in finite and infinite dimensional spaces, one may 
consult the monographs of Christensen [5] and Han et al. [13]. 

3. Sharp upper bounds on eigenvalue sums 

A symmetry of D is an orthogonal matrix U that maps D to itself (with 
the matrix acting on the left by x i— > Ux) . A group IA of orthogonal matrices 
is irreducible if every nontrivial orbit spans M. d , that is, if {Ux : U G hi} spans 
W 1 for every x £ M. d \ {0}. Equivalently, the group is irreducible if the only 
subspaces of M. d that are invariant under the action of the group are M. d and 
the zero subspace. 

Dirichlet and Neumann boundary conditions. The next theorem shows 
how a linear transformation affects the eigenvalues of a domain having irre- 
ducible symmetry group. 

Theorem 3.1. If the symmetry group of D is irreducible, then 

(Ai H h A„)| T(D) < ^[|r _:L ||| fS (Ai H h K)\ D 

for each n > 1 and each invertible linear transformation T of Mr. Equality 
holds if T is a scalar multiple of an orthogonal matrix. 
The same inequality holds for the Neumann eigenvalues. 

The theorem and its corollaries below are proved in Section 5. The two 
dimensional version of the result was proved already in our previous paper 
[20, Theorem 3.1]. In two dimensions, the irreducibility hypothesis simply 
means D has A r -fold rotational symmetry for some N > 3. Further, in two 
dimensions we found necessary and sufficient conditions for equality for the 
fundamental tone (n = 1). Here in Theorem 3.1 we state a sufficient condi- 
tion for equality, but we have not succeeded in finding necessary conditions. 

Stretching a domain by the same factor in all directions causes the eigen- 
values to scale monotonically. The next corollary obtains a monotonicity 
estimate when the domain is stretched by different factors in different direc- 
tions. Such estimates are especially interesting in the Neumann case, where 
domain monotonicity is unavailable. 

Corollary 3.2 (Stretching and monotonicity). Suppose the symmetry group 
of D is irreducible. Let T be a diagonal matrix with entries ii,...,t<f > 0, 
so that T{D) stretches D by a factor t{ in the ith coordinate direction. Then 

2 2 

(Ai + • • • + \n){T(D)) < -i + " d +td (Ai + • • • + X n )(D), n > 1. 

(Equality holds if t% = ■ ■ ■ = t^-) This result holds also for the Neumann 
eigenvalues. 
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In particular, ifti > 1 for each i, then the eigenvalue sums of the stretched 
domain are smaller than the eigenvalue sums of the original domain. 

The next corollary expresses our theorem in more geometric terms. We 
state it only for linear images of certain regular solids, for simplicity, al- 
though analogous results hold for any domain D having an irreducible group 
of symmetries, such as Archimedean solids and demihypercubes. 

Corollary 3.3 (Simplexes, parallelepipeds, and ellipsoids). Let D be a reg- 
ular simplex (respectively: hypercube, ball) in M. d . Among all simplexes 
(respectively: parallelepipeds, ellipsoids) of the form T(D), where T is an 
invertible linear transformation, the quantity 



yl+2/d 

(A 1 + --- + A n )^ 2/d 



T(D) / 



T-t(D) 



is maximal for the regular simplex (respectively: hypercube, ball), for each 
n > 1. 

The same result holds for Neumann eigenvalues. 

Other maximizers are possible too, in some cases. In particular, each 
rectangular box centered at the origin is a maximizer when n = 1 and D is 
a hypercube, as we will show in an example below. 



Understanding the geometric factors. An unusual feature of the corol- 
lary is that the eigenvalues are computed on one domain, T(D), while the 
scale-invariant geometric factor V 1+2 ^ d /I is computed on a different do- 
main, T~\D). The Example below explains why it is natural to evaluate 
the eigenvalues and the geometric factor on different domains. 

The geometric factor could be evaluated on T^ 1 (D) (without the trans- 
pose), because the transpose changes neither the volume nor the second 
moment of mass of the domain (see Lemma 5.2 and note that T _1 and its 
transpose have the same Hilbert-Schmidt norm). We have chosen to state 
the "transposed" version of the corollary because it leads to a conjecture for 
polar duals, in Section 4. 

In two dimensions, Corollary 3.3 reduces to saying that (AH hA n )^ 3 /I 

is maximal for the equilateral triangle among triangles (respectively: square 
among parallelograms, disk among ellipses), as one sees by taking d = 2, 
V = A, and using an identity connecting the moments of inertia of T(D) and 
T^ 1 (D) in two dimensions [20, Lemma 5.4]. This two dimensional version 
of Corollary 3.3 appeared in our earlier paper [20, Corollary 3.2]. 

An obvious generalization of this two dimensional estimate would be to 
maximize (Ai + • • • + \ n )V l+4: l d / 1 , but this functional is unbounded even 
on the family of rectangular boxes in three dimensions, as the next example 
shows. 
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Example: Rectangular boxes. Suppose D is the unit cube centered at the 
origin in three dimensions, and let T be a diagonal matrix with positive 
entries ti,t2,ts) see Figure 1. Then T(D) is a rectangular box with side 
lengths ti,t2,ts, volume £1*2*3 and first Dirichlet eigenvalue 

Al =7T 2 (^ 2 +t2 2 + ^ 2 )- 

And T~^(D) is a box with side lengths t^ 1 ,^ ,t7 , volume (*i*2*3) _1 arid 
second moment of mass 

(q 2 + t 2 2 + t 3 2 )/i2t 1 t 2 t 3 . 

The sum of reciprocal side length squares appears in both these last two 




i = (ti + q + t 2 3 )v/i2 i = (t- 2 + 1- 2 + tg 2 )v/i2 



Figure I. The first eigenvalue Aj = vr 2 (t 1 2 + t 2 2 + * 3 2 ) 
of T(D) is compatible with the second moment of mass of 
T~^(D), not the moment ofT(D). 



12tt z 



displayed formulas. These factors cancel to give 

v 1+2 / 3 

XiV 2/3 — 

T(D) I 

regardless of the side lengths ti,t2,ts. Thus: 

every rectangular box maximizes the first Dirichlet eigen- 
value (n = 1) in Corollary 3.3. 

Now choose t\ = t 2 = 1 and t 3 = e, so that the box T{D) is a cube that 
has been squashed in one direction. We calculate 



Ai 



T/l+4/3 



I 



12vr- 



.(2 + 6- 



.-4/3 



T(D) 



2 + e 1 



which tends to oo as e — > 0. Hence normalizing the first eigenvalue by the 
volume and moment solely of T{D) fails to yield a bounded quantity, even 
on the class of rectangular boxes. 

The preceding example helps explain why the second moment of mass 
should be evaluated on the "inverse" of the domain on which the eigenvalues 
are evaluated. Or to put it heuristically, the eigenvalue Ai is dominated by 
contributions from the "short" directions of T(D). These short directions 
contribute little to the moment of mass of T{D) but contribute a lot to the 
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moment of mass of T '(D), since they correspond to "long" directions in 
that domain . 

Robin boundary condition. 

Theorem 3.4. If the symmetry group of D is irreducible and T is an in- 

vertible linear transformation ofM. d , then 

1 



(Pl + "- + Pr0L||T-l|| Hg /Vd,T(O) ^ lW T 1 H!fs(Pl + ---+Pn)L ; 



D 



d> 

for each n > 1. Equality holds if T is a scalar multiple of an orthogonal 
matrix. 

The subscript a, D indicates both the Robin constant a and the domain 
D. Incidentally, the Robin parameter on T(D) in the theorem is multiplied 
by ||T -1 ||#s to ensure that the eigenvalues scale correctly with respect to 
T. 

A simpler and more "geometric" estimate on the Robin eigenvalues can 
be deduced by fixing the volume of the domain, as in the next corollary. 

Corollary 3.5. If the symmetry group of D is irreducible and T is an 
invertible linear transformation with | detT| = 1, then for each n > 1, 



yl+2/d 

{ P l + ---+Pn)V 2 l d 



a, T(D) I 

is maximal when T is an orthogonal matrix. 



T-t(D) 



Notice the Robin parameter in this corollary is independent of the trans- 
formation T. 

Flat torus. Our final result concerns the eigenvalues of a flat torus. Sup- 
pose T is an invertible linear transformation on M d , so that TZ rf is a lattice 
and R d /TZ d is a flat torus. Denote the eigenvalues of the Laplacian on this 
torus by Tj , so that 

= T\ < t 2 < t 3 < . . . 

and the eigenfunctions satisfy — Auj = TjUj on the torus. Some extremal 
results on those eigenvalues were obtained by Nadirashvili et al. [17, 22]. 
These eigenvalues are given explicitly as 4-7T 2 times the squares of lengths of 
vectors in the dual lattice T~^Z d , because each vector y in the dual lattice 
generates an eigenfunction e 27Tty ' x . 

Proposition 3.6. Among all flat tori R rf /TZ d , where T is an invertible 
linear transformation, for each n > 2 the normalized eigenvalue sum 

(r 2 + ---+r w ) | 2 
W 1 Whs 

is maximal when T is a scalar multiple of an orthogonal matrix (that is, 
when the torus is cubical). 
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The Hilbert-Schmidt norm in this corollary measures the "size" of the 
fundamental parallelepiped T _ t([0, of the dual lattice T~^7* d , since the 
norm is computed from the sum of the squares of the lengths of the column 
vectors of T~^. Thus the corollary can be interpreted as maximizing a ratio 
of sums of squares of lengths of vectors in the dual lattice: the numerator 
is the sum of squares of the lengths of the n shortest vectors in that lattice, 
while the denominator is the sum of squares of the lengths of d edges of the 
fundamental parallelepiped. To obtain the best possible bound, for a given 
lattice, we should choose T so that the chosen fundamental parallelepiped 
has minimal sum of squares of spanning vector lengths. 

4. Open problem — general convex domains 

We want to extend our eigenvalue bounds to general convex domains, gen- 
eralizing from the class of linear images of rotationally symmetric domains. 

The polar dual of a bounded convex domain Q in W 1 containing the origin 
is the bounded convex domain 

n° = {x £ R d : x ■ y < 1 for all y £ fi}. 

The polar dual of a ball of radius r centered at the origin is a ball of radius 
1/r. The polar dual of a square is a square rotated by 7r/4. The dual of a 
regular polygon is a rotated regular polygon. However, in three dimensions 
the polar dual of a cube centered at the origin is an octahedron. Hence the 
polar dual need not preserve the shape of a regular solid. 

If T is an invertible linear transformation then the domain T(f2) has polar 
dual T - t(f2 ). For more on polar duals, see [29, Sec. 2.8]. 

These examples suggest that taking a polar dual is similar to passing 
from T(D) to T~'(D). Therefore it seems reasonable to try maximizing the 
quantity 

yl+2/d 



{\ 1 + ... + \ n )V 2 / d 



I 



(4.1) 



This quantity is certainly bounded above, as follows. An ellipsoid E can be 
constructed inside Q such that a rescaled ellipsoid dE contains (John's 
ellipsoid [12, Chapter 4.5 Section 2]). Then (dE)° C C E°. Therefore 
quantity (4.1) can be estimated above and below, using domain monotonic- 
ity, by a constant times the same quantity for E. And (4.1) for the ellipsoid 
E is bounded above by the value for the ball, by Corollary 3.3. Thus the 
expression (4.1) is bounded above on the class of convex domains containing 
the origin. 

Problem. Find convex domains containing the origin that maximize (4.1), 
especially for the fundamental tone (n = 1). 

For Neumann eigenvalues, boundedness of the counterpart of (4.1) follows 
from the relation fj,j < Xj between the Neumann and Dirichlet eigenvalues. 
As to a maximizer, we suggest: 
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Conjecture 4.1. Suppose C K d is a bounded convex domain containing 
the origin. Then the scale-invariant Neumann eigenvalue sum 



(M2 + • • • + Vn)V 2/d 



yl+2/d 



is maximal when ft, is a ball, for each n>2. 

In two dimensions, the Conjecture claims that the product 

^421 

Q I 



(//2 H h ^n)A 



(4.2) 



is maximal for the disk. When Q is a centered ellipse, parallelogram, or 
triangle, this last product equals the simpler-looking quantity 



A 3 

(jt2 + ~'+Hn)-y 



(4.3) 



by Lemma 5.3 below. Thus for these special domains in two dimensions, 
the Conjecture is equivalent to Conjecture 4.2 in [20]. For general convex 
domains, the quantities (4.2) and (4.3) are at least comparable above and 
below (since the ratios A 2 /I on Q and f2° are comparable, using John's 
ellipsoid once more), but the relation between the two conjectures is unclear 
to us. 



5. Proofs 

The Rayleigh quotients for each boundary condition are defined as follows 

f D \Vu\ 2 dx 



Dirichlet: 
Robin: 
Neumann: 



R[u] 



j D u 2 dx 



for u G Hq(D), 



r , fn |Vu| 2 dx + £7 f an U 2 ds i 

R[u} = ^ — r — for u G H (D), 

J D u 2 dx 



R[u] 



j D \Vu\ 2 dx 
J D u 2 dx 



for it G H\D). 



The Rayleigh Poincare Principle characterizes the sum of the first n > 1 
Dirichlet eigenvalues as: 



Ai + • • • + A n 

= min{i?[ui] H + R[v n ] : v\, . 



G Hq(D) are pairwise orthogonal in L 2 (D)}. 



and similarly for sums of Neumann or Robin eigenvalues, using trial func- 
tions in H 1 instead of H]. (See [2, p. 98].) 
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Tight frame identities. Our proofs will use a Plancherel-type "tight frame" 
identity. Before stating the identity, we note that the symmetry group of 
a domain in euclidean space is compact, as one sees by a short argument 
(using that the symmetry group of the domain equals the symmetry group 
of the complement of the domain, which is a closed set). 

Let IA be the symmetry group of a domain D, so that U is compact. Write 
v for the Haar probability measure on U. If the group is finite then Haar 
measure is simply the normalized counting measure. 

Lemma 5.1. Let U be the symmetry group of a domain D. Then U is 
irreducible if and only if for all row vectors z E M. d and all matrices Y with 
d rows one has 

JjzUY\ 2 dv(U) = ±\z\ 2 \\Y\\ 2 HS . (5.1) 
In case IA is finite, the last formula says 

^£i^i» = ±i«rW™. (5.2) 

1 1 ueu 

If the matrix Y has only one column and has unit norm, then equation 
(5.2) reduces to the tight frame property (2.1) for the system of unit vectors 
{UY : U G U}. This case was proved by Vale and Waldron [26] (see also 
[27, 28]). As they remarked, the result is essentially Schur's Lemma, and we 
give such a proof below. 

Proof. Note that 

zUY\ 2 du{U) = [ {zUY)(zUY)^ du{U) 
u Ju 



z 



j ' UYY j U j dv{U)\ z ] = zMz\ (5.3) 



say, where the matrix M is clearly symmetric. Furthermore for any U £U, 
we see MU = UM by the group property of U. Let a be an eigenvalue of 
M with eigenvector w. Then 

M(Uw) = U(Mw) = U(aw) = a(Uw). 

Hence the entire orbit {Uw : U G U} consists of eigenvectors belonging to a. 
The orbit spans all of R , by the irreducibility hypothesis, and so M = a Id. 
Taking the trace yields 



ad = tvM= [ ti(UYY^U ] )dv(U) = [ tr(YYt) du(U) = \\Y \\ 2 HS . 
Ju Ju 



Solving for a and substituting M = aid. into (5.3) proves (5.1). 

"<£=" If U is reducible then for some vector y / (that is, for some 
nonzero matrix Y having d rows and 1 column) the orbit {Uy : U € U} does 
not span M. d , and so some nonzero vector z is orthogonal to every element 
of the orbit. The left side of (5.1) then equals zero while the right side does 
not, meaning (5.1) fails for this z and Y. □ 
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Proof of Theorem 3.1. The proof goes like in two dimensions [20, The- 
orem 3.1], except that the tight frame identity is more sophisticated. We 
provide a complete proof, for the reader's convenience. 

We will prove the Dirichlet case. The idea is to obtain trial functions on 
the domain T(D) by linearly transplanting eigenfunctions of D, and then to 
average with respect to the rotations of D. The Neumann proof is identical, 
except for using Neumann eigenfunctions. 

Let u±, U2, 1*3, • • • be orthonormal eigenfunctions on D corresponding to 

the Dirichlet eigenvalues Ai, A2, A3, Consider a symmetry (orthogonal 

matrix) U that fixes D. Define trial functions 

Vj = Uj o U o T _1 

on the domain E = T(D). 

The functions Vj are pairwise orthogonal, since 

Je Jd 
if j 7^ k. Thus by the Rayleigh-Poincare principle, we have 

j=l 3=1 JE V 3 ax 

For each function v = Vj we evaluate the last Rayleigh quotient as 

f E \Vv\ 2 dx _ ^ D \{Vu){x)UT~ x \ 2 dx -\detTU- l \ 
f E v 2 dx f D v? dx -\ det TU~ l \ 

= [ \(Vu)UT- l \ 2 dx, 
Jd 

where the gradient Vu is regarded as a row vector and in the last line we 
used that u = Uj is normalized in L (D). 

By averaging the preceding equality over all symmetries U in the group 
IA of symmetries of D, we find 

n . n „ 

J2\j(E)< / V / \{yuj)UT- l \ 2 dxdu(U) 
i= i Ju j=1 Jo 

iViijPUT" 1 !!^ j dx by Lemma 5.1 



j 

n 



which proves the inequality in Theorem 3.1. 

Obviously equality holds in the theorem if T is orthogonal, since the 
eigenvalues of the Laplacian are invariant under orthogonal transformations 
and the Hilbert-Schmidt norm of an orthogonal matrix equals yd. Equality 
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also holds for scalar multiples of orthogonal transformations, by a simple 
scaling argument. 

Geometric interpretation of the Hilbert Schmidt norm. To prove 
Corollary 3.3 we must express the Hilbert-Schmidt norm of T _1 in terms of 
volume and the second moment of mass. 

Lemma 5.2. If D has an irreducible symmetry group andT is an invertible 
d x d matrix, then 



In two dimensions we were able to develop a simpler formula [20] solely in 
terms of the domain T(D), due to a relation between the Hilbert-Schmidt 
norms of T and its inverse (||T||#s = | detTlUT -1 !!^). No such relation 
seems to hold in higher dimensions. 

Proof of Lemma 5.2. The centroid of D lies at the origin, as a consequence 
of the irreducibility of the symmetry group. Hence the centroid of T~\D) 
also lies at the origin, by a linear change of variable. 

The moment matrix of D is defined to be M(D) = f D xx^ dx, where x 
is a column vector. We claim M(D) is a scalar multiple of the identity. 
For let U be a symmetry of D. The invar iance of D under U implies that 
M(D) = UM(D)U\ so that M is a scalar multiple of the identity by arguing 
like in the proof of Lemma 5.1 (that is, Schur's Lemma again). 

Since the diagonal entries of M(D) are equal, they must be 1/d times the 
trace of the moment matrix, that is, 1/d times the second moment of mass 
(using here that the centroid of D lies at the origin). Thus 




M(D) = -1(D) Id. 



(5.4) 



The moment of inertia of T ^(D) can now be computed as 



I(T^D) = tr M(T~^D) 

= tr (r _t M(D)r _1 | detT" 1 !) 

= -1(D) (trT^T^IdetT- 1 ! by (5.4) 
= i/(D)||T- 1 ||^|detT- 1 |. 



Rearranging, 



d" UHb | det T- 1 ! 1(D) 



V-iii2 i J(r-tg) 



(5.5) 



The lemma now follows easily, using the formula 
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□ 

Next we divert from our main argument to observe that certain plane 
domains deviate from roundness to the same extent as their polar duals, 
when the deviation is measured by the moment of inertia. 

Lemma 5.3. If Q is an ellipse, parallelogram, or triangle, and the centroid 
of fi lies at the origin, then 



I 



Proof of Lemma 5.3. Write D for a disk, square, or equilateral triangle, 
with centroid at the origin, in dimension d = 2. Then by the examples 
in Section 4, we know the polar dual D° has the same shape as D, up to 
rotation and scaling. Hence the scale invariant ratio A 2 /I takes the same 
value for D° as for D. 

The domain fi is a linear image of D, meaning = T(D) for some 
invertible linear transformation T. We compute 

^(t(d)) = ^r'lmf b y ( 5 - 5 ) with T re P laced b y T ~ ] 

1 W 1 Whs i \ IJ ) 

_ 2| detTl" 1 A{D) 2 

~ \\T- l f HS HP) 

because HT -1 !!^ = detT| 2 in two dimensions 



2|detT" 1 | A(D°] 
WT-'Whs !(D° 

y(rt(D°)) 



2 



since D and D° have the same shape, 



by (5.5), which proves the lemma because T~\D°) = T(D)° = □ 

Proof of Corollary 3.3. For the ball, irreducibility of the symmetry group 
is obvious since every rotation is a symmetry of the ball. If D is a regular 
simplex or a hypercube, then the symmetry group is known to be irre- 
ducible; for example, see [6, Chapter 6], where irreducibility is shown for 
the symmetry groups of the Platonic solids and higher dimensional regular 
polytopes. 

Corollary 3.3 therefore follows from Theorem 3.1 and Lemma 5.2. 

Proof of Theorem 3.4. The Robin case in higher dimensions is signifi- 
cantly more involved than in the planar case, because to handle the bound- 
ary term in the Robin Rayleigh quotient we must study the Jacobian of T on 
hyperplanes tangent to the boundary of D. (In two dimensions the bound- 
ary has only one tangent direction, and the boundary Jacobian is easily 
understood.) 

We may assume the linear transformation T is diagonal with positive en- 
tries t\, . . . ,td, since the general case may be reduced to this situation by the 
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singular value decomposition of T and the invariance of Robin eigenvalues 
under rotations of the domain. 

We will study the action of T on a hyperplane by examining its action 
on the normal vector, motivated by the fact that the base area of a par- 
allelepiped can be determined from the volume and the altitude. Consider 
vectors wx,W2, . . . , Wd-i and w. Define W = [wi . . . Wd-i] to be the matrix 
with iih column Wi, and let 



S[W, w] = det[u>i . . . Wd-i w}. 

Write e for the column vector whose ith component is e^, the ith standard 
unit column vector. Define 

S[W] = S[W, e] = det[wi . . . w d -x e], 

so that S[W] is a vector perpendicular to each vector u>j. Notice I.SfW 7 ] • 
w\ = \S[W,w]\ gives the volume of the parallelepiped spanned by vectors 
wi, . . . ,Wd-i,w. By choosing w to be a unit vector orthogonal to each Wi, 
we deduce that |iS[W]| equals the (d— l)-dimensional volume of the (d— 1)- 
dimensional parallelepiped spanned by w\, . . . , Wd-i- 
One can check from the definition that 

S[TW] = 5[^,T- 1 e](detT). (5.6) 

Further, for any orthogonal matrix U one has 

US[W] = S[UW] (5.7) 

by the geometric interpretation of S[W] as a vector perpendicular to each 
column of W and with magnitude equalling the (d — l)-volume of the par- 
allelepiped spanned by wi, . . . , w d -i- 

Now we can prove the theorem. The proof goes like for the Dirichlet 
and Neumann cases in the proof of Theorem 3.1, except that for the Robin 
eigenvalues we must take account also of a boundary integral in the Rayleigh 
quotient. Following the notation of that proof, the boundary term is 

J dE v 2 dS(x) _ j dE u(UT^xfdS{x) 



J E v 2 dx J E u(UT l x) 2 dx 

_ f 9D u(x) 2 \S[TU-'W(x)]\dS(x) 
f D u(x) 2 dx ■ | det T\ 

by a change of variable, where W(x) = [wi...w d -i] is a matrix whose 
columns form an orthonormal basis for the tangent space of D at s. Hence 



u(x) 2 \S[U- 1 W(x),T- 1 e}\dS(x), (5.8) 

dD 



J E v 2 dx 

by calling on (5.6) and using the normalization of u in L 2 (D). 
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The Jacobian factor in (5.8) satisfies 

d 

\S[U- l W,T- l e}\ 2 = 



s[u- l w\ 



i=i 



8=1 



t 2 



by (5.7). Integrating over all symmetries U and then applying the tight 
frame property from Lemma 5.1 shows that 

f ISiU-^T^e}] 2 dv(U) = / {ei-U^SiW}) 2 du(U) 

Ju i^i t i 



i=i 1 

W L Whs 

d ' 



\S[W}\ 



since orthonormality of the vectors wx, . . . , Wd-i guarantees that |i?[W]| = 1. 
Hence by taking the square root and using Cauchy-Schwarz, 

/ \S[U- 1 W,T- 1 e]\dv(U) < J!ZLi^£. 
Ju yd 

Applying this last estimate to (5.8) gives that 



/ 

Ju 



J dE v 2 dS(x) 



du{U) < 



\HS 



u(x) 2 dS(x). 



IU f E <' 2 'I-'' \'d Jill) 

The adaptation of the proof of Theorem 3.1 to the Robin situation can 
now be completed without difficulty. 

Proof of Corollary 3.5. The inequality between quadratic and geometric 
means implies that 



\HS 



Vd 



— (singular value of T 1 ) 2 



i=i 



> 



\ ]^[( s i n g iuar value of T 1 ) ? : = \J | det T 1 1 
\ i=i 



by our hypothesis that |detT| = 1. By inserting this last inequality into 
Theorem 3.4 and invoking the monotonicity of the Robin Rayleigh quotient 
(and hence eigenvalues) with respect to the Robin parameter, we find 



1 



[ftl + V pn) 

Now Lemma 5.2 completes the proof. 



< -\\T 



Pr, 
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Proof of Proposition 3.6. Write Vi, ... ,v n for a collection of n shortest 
vectors in the cubical lattice Z d (so that v\ = 0). Then since T~*vi, . . . , T^ Jf v 1 
are distinct vectors belonging to the dual lattice T~^Z d of TZ rf , we have 

(ri + • • • + r n )| Rd/TZd < 4tt 2 (\T^ Vl \ 2 + ■■■ + \T^v n \ 2 ^ . 

Let U be any symmetry of the hypercube [—1/2, l/2] d , so that U maps Z, d 
to itself. Then another collection of n shortest vectors in Z rf is Uv\, . . . , Uv n . 
Repeating the above estimate with this new collection, we find 

(ri + • • • + r n )\ Rd/TZd < in 2 (\T^Uv x \ 2 + ■■■ + \T^Uv n f} . 

Averaging this inequality over all symmetries U of the hypercube implies 
(by the tight frame identity in Lemma 5.1) that 

(ri + • • • + T n )\ Rd/TZd < Air 2 (M 2 + • • • + \v n \ 2 ) i! 

= (n H h T"n)| Rd / Z d |, T ,_1||2 • 
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